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[TECHNICAL NOTE k 230 


PRANDTL-MEYER EXPANSION OF CHEMICALLY REACTING GASES 
IN LOCAL CHEMICAL AND THERMODYNAMIC EQUILIBRIUM 
By Steve P. Helms 


SUMMARY 


It is found that Prandtl-Meyer flow, in which chemical reactions * 
are occurring and ar e in e quilibrium, can he simply and exactly calcu- 
lated. The property of" air wmch governs the flow is found to he a 
quantity which depends only on the ratio of enthalpy to the square of 
the speed of sound; the analogous quantity for an inert gas depends only 
on the ratio of specific heats. The maximum angle through which the 
flow may turn is generally larger when chemical reactions are occurring - - 
than it is in nonreacting air. A numerical example shows that the pres- 
sure variation with angle, as well as temperature and Mach number varia- 
tions, may he considerably affected hy the presence of the chemical 
reaction. 


INTRODUCTION 


At the high temperatures encountered in hypersonic flight the air 
may no longer he regarded as an inert gas. It does not have a constant 
ratio of specific heats, 7, nor does it generally obey the simple equa- 
tion of state, p/p = RT. These thermodynamic features reflect the fact 
that at high temperatures molecular vibrations are excited and chemical 
reactions are taking place in air. Because of this, any flow solutions 
depending on constancy of 7 and the perfect gas law are not valid. 

One elementary supersonic flow solution is the Prandtl-Meyer expan- 
sion around a comer. In this paper the theory of the Prandtl-Meyer 
expansion is extended to include high temperature flow in chemical and 
thermal equilibrium. When only molecular 1 vibrations are active, and no 
chemical reactions occur, one can still use the usual flew equations in 
terms of 7, if the appropriate function of temperature is inserted for 
7 (see ref. l) . However, when chemical reactions are occurring, then 7 
is no longer a useful concept in the Prandtl-Meyer flow. 

Instead of working with 7, we shall employ a quantity tj to des- 
cribe the thermodynamics of the gas, because it is q and not 7 which 
enters into the flow equations at high temperatures. When no reactions 
are taking place, q reduces to (7+l)/(7~l). In the present analy sis, 
it will be shown that by introducing also an auxiliary variable lfr, the 
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Prandtl -Meyer expansion can be simply and exactly calculatec for equi- 
librium flow without employing any iterative procedures or extensive 
numerical integrations. However, a table orTa Mollier diagram for the 
thermodynamic properties of air at high temperatures is required. At 
high temperatures, the table of reference 2 or the Mollier chart of 
reference 3 is suitable. At temperatures below 3000° K, reference ^ is 
useful. 


SYMBOLS 


a speed of sound 

c a constant of -motion with dimensions of velocity, defined by 
equation ( 3 ) 

h enthalpy per unit mass 

jd Planck’s constant divided by 2jt 

k Boltzmann constant 

M Mach number 

p pressure 

r radial coordinate 

R gas constant for air 

S entropy 

T absolute temperature 

v speed, *Jv r 2 + Vg 2 

v r component of velocity along the radius vector 

v _ component of velocity perpendicular to the radius vector 

9 

c p 

7 ratio of specific heats , 

2h 

n parameter defined by 1 + — 

or 

G angular coordinate 

v angle through which the flow has turned 

v m theoretical maximum value of v 

p density 
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i|r an angle defined by equation (7) 
an integration variable 

Subscript 

o value along the radial line on which the Prandtl -Meyer expansion 
begins 


ANALYSIS 


In the usual derivation of Prandtl -Meyer flow (ref. 5 ) one seeks a 
solution for supersonic flow for which the pressure, density, and velocity 
are constant along radial lines emanating from the corner around which 
the gas is expanding. Such a solution can exist only when no character- 
istic length enters into the problem. Thus when chemical relaxation 
occurs in the air, one cannot have such a solution. However, when the 
air is in equilibrium everywhere, there is no length and one may expect 
to find such a solution. More generally, if each of the various reac- 
tions and internal degrees of freedom of the molecules are either frozen 
or in chemical equilibrium, a Prandtl -Meyer type solution is expected to 
exist. We seek a formal solution of the Prandtl-Meyer type which is 
valid for any of these isentropic flows. 


If we define 



the variation of pressure with density at 


constant entropy, we obtain from the requirement of continuity and from 
momentum conservation the usual equations: 



CD 


v e = a 


( 2 ) 


* 


* 


Equation (2) requires that the flow be supersonic. The energy equation 
is 


1 

2 


h + 


2 2 
V„ + V 

9 r 


Combining equations (l), (2), and (3) gives the differential equation 
for v r : 


( 3 ) 
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(a/) “ 1 + (2h/ a 2 ) (° 2 " v r 2) 

Let 0 o be the angle at which the Prandtl -Meyer expansion begins. Then 
in view of equation (2), the initial radial velocity is (see sketch): 


v r(®o) ~ */Mq 2 ” 1 


( 5 ) 


8-0 



It is seen from equations (l) and (!<■) that, for given initial condi- 
tions, the flow is completely specified if the quantity rj is known, 
where we define 


il = 


l + % 
a 


(6) 


For air at moderate temperatures, the. ratio of specific heats is a 
constant, and t} = (7+1) / (7-1) = 6. For a diatomic gas in which the 
molecular vibrations are in equilibrium but no reaction is occurring, T| 
is a function of temperature (but not of pressure) and varies from 6 at 
moderate temperatures to 8 at temperatures that are large compared to 


* 


* 
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]&o/k, where to is the characteristic frequency of molecular vibrations. 
In general, when chemical reaction is also occurring, the quantity tj 
is a function of local temperature and pressure. 

In order to integrate equation (4) in this general case, it is 
convenient to introduce a variable Tjf(O) defined by 


dg, _ /- 

dijr " 

or > 

rt 

0 - 0 O - / -s/n af 

J 

with 

t 0 = 8ln ®o) = sin " 1 (^V M o 2 “ l) 

Then equations (4) and (5) become 



(7) 


(8a) 


v r(to) - ®o “ 1 


(8b) 


The solution of equation (8) is 


v r - c sin t , t > i(r 0 


(9) 


The other component of velocity is then in view of equations (l) and (2) : 


— 1/2 

= a = ct} cos ijr 


( 10 ) 


It is seen from equations ( 9 ) and (10) that when the angle i(r Is 
it/2, the flow is purely radial and the Mach number becomes infinite. 
From this it follows that jt/2 is the limiting value of ljr. 

To compute the flow field, we notice from equation (10) that ijr 
can be expressed as a function of the state of the gas: 
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♦ = c ° s_ 1 (!M) (n) 

Now with equations (2), (6), ( 9 ), and. (ll) , we can compute the flow field 
as a function of if. When chemical reaction is important the simplest 
procedure is to use a Mollier chart for air (ref. 3 ) or a set of tables 
(ref. 2), and regard h as the independent variable. Knowing the initial 
(\|r = + 0 ) state of the gas, and knowing that it will follow an i sent rope, 
one takes successive increments of h, looks up corresponding values of 
the sound speed a, readily evaluates the right-hand side of equation (ll), 
and thus finds the value of \|r corresponding to each h. Any of the 
other thermodynamic properties of the air corresponding to a particular 
can then also be read off from the chart or table. At the lower tem- 
peratures where the molecular vibrations are active, but where the chemi- 
cal reaction may be neglected, a somewhat simpler procedure is suggested 
by the fact that tj and a are in this region functions of temperature 
only: Choose temperature as the independent variable (an is done in the 

tables of ref. 3)> and look up h and a in the tables. Then equation ( 11) 
yields the value of \|f . 

In any case, once i|f is known, one obtains the flow field from 
equations ( 9 ) and (10), and the local Mach number is given by 


M = + (v r /vg) 2 = J 1 + q tan 2 ^ 


( 12 ) 


It remains to obtain the physical angle 9 corresponding to the parameter 
if by performing the integration indicated in equation (7 ) . At tempera- 
tures where molecular vibrations and chemica l reaction ca n be neglected, 
this integral (for air) is simply Q - 0 o = n/ (7+ 1)/ ( 7 - 1 ) (if - ijr 0 ). Even 
at high temperatures varies only slowly (as fig. 1 illustrates), and 

numerical evaluation of the integral ( 7 ) can be done quickly because one 
may take large steps, treating as constant during each step, 

This completes the solution of the Prandtl -Meyer problem; the angle 
the flow has turned through is given by: 


v = 0 


cos 



( 13 ) 


(For an alternative expression see ref. 3*) 

For some applications it may be desirable to have an explicit analytic 
expression for the density an a function of ijf. Such an expression can be 
obtained by integrating the continuity equation, (dpvg/dS) + pv r = 0 (see 
ref. 5). The integration is performed with the help of equations ( 7 ), ( 9 ), 
and (10). It yields: 
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p(t) 

p(*o) 


Ot 

/ , \n - 1 / m/2 / (v T >) taa d * 

/ cos Tjr yto YJLY e ^o 
vpos to/ \Jlo/ 


(14) 


The usual "constant 7 " solution is obtained when t) = tj q = ( 7 + 1 )/ ( 7 - 1 ) . 
A numerical table of h and a permits fitting (tj q - tj ) / ti 0 , the local 
deviation from constancy, by a simple sine series with integer multiples 
of t ,- t 0 as the argument. Ihe integral in equation (l4) is then of an 
elementary type. In the special case where Mq is unity, the integral 
in equation (l4) may be put equal to zero for small angles ijr. 


NUMERICAL EXAMPLE AND GENERAL QUALITATIVE FEATURES 


A numerical example of Prandtl -Meyer flow is presented graphically 
in figures 1 to 5* A set of initial conditions were chosen at which the 
difference between frozen and equilibrium flow is marked, and at the same 
time is in the region which is of interest for high-altitude flight. It 
corresponds to the stagnation conditions at an altitude of about 130,000 
feet with a flight velocity of 19 j 000 feet per second. The frozen flow 
was taken as a flow in which only the random translations and molecular 
rotations keep up with equilibrium, but the chemical reaction and the 
energy of the molecular vibrations are frozen at the values they have at 
the point where the Prandtl -Meyer expansion begins. The initial thermo- 
dynamic state (pressure , temperature, density) at a Mach number of unity 
is taken to be the same in the frozen and equilibrium case . This requires 
that the initial sound speed and the flow velocity are different in the 


frozen case than in the equilibrium case. 


The sound speed, 



is 


different in the two cases because the variation of pressure with density 
is different, even though the initial thermodynamic states are identical. 
This difference in sound speed results in a significant difference in the 
initial value of T] (see fig. 1 ) . 


Several qualitative features of the relation between equilibrium 
flow (large q) and frozen flew (small q) can be seen from the present 
analysis and from the figures. Because of the simple geometry, we 
restrict ourselves in the following comments to the important type of 
flow in which the expansion begins with sonic velocities, even though 
al 1 of the foregoing equations are valid for any initial Mach number not 
less than unity. 


The Maximum Turning Angle 

* 


Hie maximum angle through which the flow may turn corresponds to 
i|r = 3 t /2 and is given by 


4 
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When rj is constant, equation (15) reduces to 


V m = Wn - 1) § 


(15a) 


Even when tj is not constant, equation (15a) is evidently approximately 
equal to equation (15), provided a suitably averaged value of >Jri is used. 
For frozen undissociated flow t ] - 6 and v m = 129. 3°. However, for flow 
with chemical reaction (see fig. l), the averaged value of rj may be 
larger than 9 and consequently the maximum angle v m is larger than l80°. 
This shows that at high temperatures and pressures, equilibrium flow might 
have to make more than a l80° turn before the theoretical limit is reached. 

The Mach number, pressure, density, and temperature a.l 1 take on their 
limiting values (M becomes infinite, T, p, and p become zero, for the 
ideal gas model) 1 at v = v m . This fact is reflected in all the figures 2 
to 5 in that for the larger turning angles the small q values of p, T, 
p, and M have changed more from their initial value, than the large T] 
curves. 


Density 


An interesting feature of the density curve (fig. 2) is that the 
equilibrium value is not everywhere higher than the frozen value of 
density, but its initial slope is steeper and the two curves then must 
crosB, because from the above comment about vm we know that the density 
for the frozen case has to reach zero much sooner than the equilibrium 
curve. It w ill usually be true that for the equilibrium flow the density 
decreases initially more rapidly than for the frozen flow. This can be 
proved rigorously from equation (l4) , if the equilibrium flow is identi- 


fied with a certain value of 


and a nonpositive value of 




and the frozen flow is identified with a relatively smaller and a 
constant ri 0 (see appendix). 


1 The point here is not so much that T = 0 and p = 0 at v m j in fact 
the ideal gas model breaks down completely at temperatures where the air 
begins to liquefy. For example, near T = 0 the specific heats are zero, 
and do not have the constant values associated with ideal gases at room 
temperature. Hie point is, more precisely, that the v corresponding to 
the state of the air in which the ideal gas laws break down is much larger 
in the large T) case than in the small T] flow. In the text the simpler 
statement is preferred; however, all of the equations in this paper are 
generally valid and do not depend on any specific gas model. 
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Pressure 


* First of all one sees from figure 3 that unlike the density curve 

the equilibrium values for pressure are always higher than the frozen 
ones. Moreover the equilibrium and frozen pressure curves are well sepa- 
rated, much more so than the density curves. For example at v = 70° 
the equilibrium pressure is nine times as large as the frozen pressure. 
The corresponding ratio of the densities is 2. 

The point of interest is that in a Prandtl -Meyer expansion the pres- 
sure is quite sensitive to the chemical reaction, much more so than the 
density. This is Just the opposite of what occurs in the compression by 
a normal shock, in which the pressure ratio is very insensitive to any 
chemical reaction, but the density ratio is not. 


Mach Humber and Turning Angle 


The local Mach number variation and its strong dependence on T] in 
the example are shown in figure 4. The more general qualitative variation 
of Mach number and turning angle can be derived from equations (12), (13), 
and (7). It is seen from these equations that for a given angle 0, the 
Mach number and v depend on the local value of T], as well as on its 
value for smaller 0, but they depend on no other parameters. The differ- 
ence between frozen and equilibrium flow may be analyzed as arising (a) 
from the fact that t) is usually larger in the equilibrium case and (b) 
that it is not constant in the equilibrium case. The effect (a) is seen 
from equation (12), and may be stated in words: The larger the (constant) 

value of q the smaller the value of the local Mach number at any fixed 
0, and the larger the corresponding turning angle v(0). It follows that 
the larger the value of t], the smaller the Mach number at any given 
turning angle v. As to the effect (b) due to the variation of t]( 0) 
with 0 , it may be stated in words as follows: If tj( 0) is a monotoni- 

cally decreasing function, then the actual Mach number at a point 0 X 
for equilibrium flow is smaller and v( 0 1 ) is larger than it would be 
for a flow which is characterized by the constant value of q = tj(0 1 ) . 

It follows that for any fixed v, the Mach number is smaller than for 
the flow characterized by the constant rj = t)( 0 -l) . 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Dec. 2, 1957 
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APPENDIX 


THE DENSITY AT VERY SMALL ANGLES 


We wish to consider two flows initially at M = 1 and at the same 
density (which we choose as our unit of density), hut with different 
initial values of T], t| 0 and T} 0 . Also generally the derivatives 


£? C3 
*1 - 


d log i). 


d0 


$ 


Bo*” 36 * 


d log iy 


d0 


a 


2 ) will "be different. We ask the 
9o 


question: Under what conditions will the density of flow 1 be 

smaller t han that of flow 2 for small angles? 


Expansion of equation (14) for small angles 9 yields: 


P 


1 





~ g 2 ) + 


@2 /Tlo 2 “ ^O 

2 V *0^ 



From this expression it is seen that the condition that p be less than 
P 2 for sufficiently small 9 is either 

or 

Si = Sg j n Qx ^ 

In the example given here, where p x represents equilibrium flow and p 2 
frozen flow, one has g x < g 2 = 0, as well as ti Qi > 
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Flow deflection, v 

Figure 2.- Comparison of density for frozen and equilibrium Prandtl -Meyer 
flow; Mo = 1.00, T 0 = 6l4o° K, p Q = 1.2 atmospheres. _ 



NA.CA ra 4230 


15 



0 ° 10 ° 20 ° 40 ° 60 ° 80 ° 

Flow deflection, v 

Figure 3.- Comparison of pressure for frozen and equilibrium Prandtl -Meyer 
flow; Mq = 1 . 00 , T 0 = 6l4o° K, p Q = 1.2 atmospheres. 













